Quantum Walks on a Random Environment 
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Quantum walks are considered in a one- dimensional random medium characterized by static or 
dynamic disorder. Quantum interference for static disorder can lead to Anderson localization which 
completely hinders the quantum walk and it is contrasted with the decoherence effect of dynamic 
disorder having strength VK, where a quantum to classical crossover at time tc oc transforms 
the quantum walk into an ordinary random walk with diffusive spreading. We demonstrate these 
localization and decoherence phenomena in quantum carpets of the observed time evolution and 
examine in detail a dimer lattice which corresponds to a single qubit subject to randomness. 

PACS numbers: 03.67.Lx,72.15.Rn,03.65.Yz 



I. INTRODUCTION 



During last decade quantum algorithms were proposed, 
such as Grover's search[l] and Shor's factorization [2], 
which can in principle perform certain computational 
tasks quantum-mechanically, much more efficiently than 
their classical counterparts. The related idea of quantum 
walks was also introduced^, 0, [HI, [s.] which generalize the 
classical random walks widely used in various computa- 
tions as the basis of classical algorithms. The quantum 
walks are similar to classical random walks but with a 
"quantum coin" operation which replaces the coin- flip 
randomness in between each moving step on a lattice or 
a graph. The state of the quantum coin which uniquely 
determines the subsequent movement direction can also 
exist in quantum superpositions, something impossible 
in the classical domain where the coin has a specific out- 
come. In analogy with classical random walks the quan- 
tum walks are expected to be useful for designing quan- 
tum algorithms. For example, Grover's algorithm can 
be combined with quantum walks in a quantum algo- 
rithm for "glued trees" which provides even an exponen- 
tial speed up over classical methods [7]. 

The main advantage of quantum walks is a highly im- 
proved behavior over their classical counterparts since 
quantum wave propagation is superior than classical dif- 
fusion. For example, the ballistic mean square variance 
cr^(t) oc in the quantum case can be compared to the 
linear spread law cr^(t) oc t of classical diffusion. This 
quadratic speed-up is a general feature of quantum search 
algorithms |l] and is also familiar from standard quan- 
tum evolution of tight-binding electron waves on a pe- 
riodic lattice [8]. In quantum walks the classical proba- 
bilities P(x, t) are replaced by complex probability am- 
plitudes computed from the unitary dynamics of 
the Schrodinger's equation. The corresponding proba- 
bility amplitudes are determined by summing up over 
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all possible paths of propagation. Furthermore, to de- 
scribe wave propagation in lattices or graphs one does 
not need a " quantum coin" [9] and a related continuous- 
time versions of quantum walks have been introduced[T^. 
The discrete and continuous-time quantum walks have 
recently been related to Dirac and Schrodinger's equa- 
tions, respectively [11]. 

We consider a continuous-time quantum walk via the 
equivalent problem of a quantum particle initially local- 
ized in one-dimensional lattice in the presence of static 
or dynamic disorder. For a tight-binding electron at an 
integer lattice site labeled by x in one dimension at time t 
the wave function ^(x, obeys the linear wave equation 

ih = i) 

id-^ix, t)/dt = e{x, t)'^{x, t) + -^{x -l,t) + *(a; + 1, t), 

(1) 

with e(x,t) an x-dependent random variable for static 
disorder which is also t-dependent for dynamic disorder, 
where lengths are measured in units of the lattice spac- 
ing and energies or inverse times in units of the hop- 
ping integral. In order to study quantum walks via Eq. 
(1) we choose the initial condition of a particle at the 
origin x = with ^(x,t = 0) = ^3^,0 and characterize 
the quantum motion by the second moment for its posi- 
tion cr2(t) = where P(x,t) = |^(x,t)p 
is the probability density. In the absence of disorder 
(e(x,t) = 0) the amplitude is given by the Bessel func- 
tion so that P(x,t) = (Ja^(2t))^, where the order x of 
the Bessel function measures the distance travelled from 
the origin while its argument is proportional to time t. 
The evolution of the wave-packet at time t was shown [H 
to display two sharp ballistic fronts at x = ±2t. From 
properties of the Bessel functions inside the spatial re- 
gion [— 2t,2t] the probability density P(x,t) is an oscil- 
lating function multiplied by 1/t while outside this re- 
gion, denoted by the two ballistic peaks, P(x, t) decays 
exponentially. The ballistic mean-square displacement is 

In the absence of disorder the quantum evolution of 
wave-packets which occurs via the evolution operator 
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ex.p{—iHt) with Hamiltonian H is very different from 
classical diffusion where any initial state converges to a 
Gaussian steady state. An initially squeezed J-function 
spatial wave-packet has reduced spatial uncertainty and 
behaves like a quantum particle consisting of all the 
eigenstates of H. Alternatively, a spatially uniformly dis- 
tributed initial wave-packet ^(x,0) = l/\/]V, for every 
X in an N-site chain, consists of few eigenstates of H 
near the lower band edge only. Since the latter choice 
emphasizes states from the band edge the semi-classical 
asymptotics is relevant [12[ 0, [3]. The wave-packets 
which initially consist of many eigenstates can be related 
via the uncertainty principle to ultrashort laser pulses of 
femtosecond duration. The evolution of such coherent su- 
perpositions of quantum states is realized in the physics 
of trapped atoms in optical lattices [is], trapped ionspjj, 
etc. 

Classical random walks in perfect one-dimensional lat- 
tices are defined by the probabilities Px = Qx = 1/2 
which determine the random left and right motion from 
site X. This externally induced randomization leads to 
classical diffusion for long t with a Gaussian P(x, t) and 
cr^(t) = 2Dt^ D the diffusion coefficient. For random 
walks in random media the probabilities Px^Qx = ^ — Px 
become random variables themselves, for example, they 
could be chosen from a flat distribution within (0,1). 
This is the so-called random random or Sinai walk[17l.llq 
which leads to ultra slow classical evolution cr^(t) oc Int^, 
very close to a complete cease. 

The problem addressed in this paper concerns the fate 
of quantum walks in a random environment, with both 
static and dynamic disorder. To answer this question 
for static disorder we shall combine previous knowledge 
from the field of wave propagation in the presence of ran- 
domness where the quantum phenomenon of Anderson 
localization [191 takes place (for its consequences for quan- 
tum walks see [20]). We shall show that static disorder is 
responsible for exponentially suppressed quantum evo- 
lution with variance cr^(t) reaching a time-independent 
limit for long t, depending on the strength of static dis- 
order and space dimensionality. Surprisingly, classical 
random walks for static disorder are still propagating, 
although ultraslowlvfiri Other generalizations of 

discrete-time quantum walks in aperiodic or fractal me- 
dia by using biased quantum coins have given slower 
(sub-ballistic) quantum evolution [21]. For dynamic dis- 
order by coupling the quantu ni sy stem to a random en- 
vironment decoherence occurs [22] and quantum physics 
becomes classical so that a quantum walk is still propa- 
gating but only diffusively. 

The main reason for examining the robustness of quan- 
tum walks in the presence of noise is because disorder 
is unavoidable in most quantum systems. Static disor- 
der also appears for electrons in lattices with permanent 
modifications due to impurities. The dynamic disorder 



in this case could be driven by time-dependent vibrations 
of the lattice atoms which have an impact on the elec- 
tronic site-energies. Apart from describing such electron- 
phonon interactions, dynamic disorder also addresses the 
presence of time-dependent noise in the memory qubits of 
quantum computers. In this paper we demonstrate that 
although static disorder hinders the motion of quantum 
walks due to negative quantum interference from Ander- 
son localization via multiple scattering from impurities, 
instead, for a dynamically random environment the time- 
dependent disorder acts as a decoherence mechanism at 
a crossover time tc which randomizes the quantum walks 
and turns the quantum motion into classical. The dis- 
crete finite space chosen in the simulations could be also 
used on a finite computer. This general scheme for dis- 
cretizing space is not only suggested by solid state appli- 
cations but it is, somehow, related to the discreteness of 
the quantum information itself. 

The paper is organized as follows: In Sec. I we intro- 
duced the reader to the subject of quantum communica- 
tion by setting the aims of our quantum walk simulation 
in random media. In Sec. II we briefly review the prop- 
erties of quantum walks in ballistic and disordered one- 
dimensional media by showing quantum carpets which 
demonstrate the difference between static and dynamic 
disorder. Static disorder is shown to be responsible for 
negative quantum interference of Anderson localization 
which stops completely the quantum walk while dynamic 
disorder permits only diffusive evolution of classical ran- 
dom walks. In Sec. Ill we display the quantum to clas- 
sical crossover for dynamic disorder and consider a qubit 
subject to dynamic disorder. Finally, in Sec. IV we sum- 
marize our main conclusions. 



II. QUANTUM CARPETS 

We have created space-time x — t structures for the 
probability density P{x,t) on a one-dimensional finite 
A^-site orthonormal lattice space without disorder, with 
t-independent static disorder and also rapidly varying dy- 
namic disorder. The white color in the figures denotes 
high probability density and the darker colors lower val- 
ues. In Fig. 1 a state is initially released in the mid- 
dle of the chain and the probability density P(x,t) is 
obtained by solving Eq. (1). In Fig. 2 the same is 
done for a spatially uniform initial state. For static dis- 
order Eq. (1) could be alternatively solved by consider- 
ing the time-evolution of a state vector |^(0)) expressing 
the probability density via the stationary eigen-solutions 
H\j) = Ej\j) with space-time wave function 

N 

*(x,t) = ^e-^^^*^,(x)(i|*(0)), (2) 
i=i 
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with the amphtude on site x denoted by ^|Jj{x) = {x\j). 
For time-dependent disorder Eq. (1) was solved via a 
fourth-order Runge-Kutta algorithm. 

In Fig. 1 we present our results for the initial choice of 
a ^-function in the middle of the chain with |^(0)) = |0) 
and in Fig. 2 for a uniform initial state |^(0)) = 
(I/a/TV) ^x^i 1^) computed on finite chains where the 
wave-packet scatters from their hard ends. In the rest 
of Figs. 3-8 our results are obtained for a self-expanding 
chain with a ^-function initial choice to make sure that 
the wave-packet does not reach the boundaries. This al- 
lows to study the quantum to classical crossover by com- 
puting the mean-square- variance and the autocorrelation 
function vs. time. 




ballistic motion 



In this case, obviously, quantum walks perform at 
their best. The ballistic description is valid for solid 
state systems in the absence of disorder which refers to 
the motion of a point particle in an N-site chain with 

2cos(-^^^) and 



in Eq. (1) which gives Ej 



^sin(^), j = 1,2,...,7V. The corre- 



tpjix) = ^ , 
sponding space-time pictures are shown in Figs. 1(a) 
and 2(a). Quantum revivals can be seen where the par- 
ticle returns to its initial position and reconstructs like 
a classical particle which moves with constant velocity 
refiecting at the boundaries of the chain[14]. We observe 
that the quantum revivals of Fig. 1(a) and Fig. 2(a) 
do not repeat indefinitely but become less and less ac- 
curately as time progresses. This is due to effects from 
boundary scattering which become more prominent for 
broad wave-packets in the right hand-side of Fig. 2(a) 
where noisy evolution is established. The obtained frac- 
tal pattern is a result of peculiar quantum interference 
effects due to scattering from the hard walls at the ends 
of the chain 0. 



static disorder 

It can have dramatic consequences for quantum walks, 
particularly in low dimensions, since for static disorder 
they can stop completely due to Anderson localization. 
From Figs. 1(b) and 2(b) we can see how strong static 
disorder with e{x) chosen from a uniform probability 
distribution within [— VK/2, +W/2] causes destructive in- 
terference with Anderson localization. This crossover 
from ballistic motion to localization has dramatic con- 
sequences for quantum walks in one-dimension in the 
presence of static disorder. In higher dimensions the 
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FIG. 1: Quantum carpets which show the probability density 
P{x,t), for X in the vertical axis and t in the horizontal axis, 
are generated by an initial ^-type spatial state for t — 
^(x,0) = 6x,o in the middle (left of the figure) with chain 
length N — 101. (a) The ballistic case for the absence of 
disorder where perfect quantum revivals can be clearly seen, 
(b) For static disorder of strength VK = 1.5 quantum interfer- 
ence causes Anderson localization which stops the quantum 
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(a) 




FIG. 2: Quantum carpets showing the probabihty density 
P{x,t) for a spatially uniform initial state ^(x,0) = 
on a chain of length N — 101. (a) For the ballistic case in 
the absence of disorder one can see recurrences for short times 
which disappear for longer times where quantum interference 
effects become apparent, (b) For static disorder of strength 
W — b Anderson localization occurs with localization length 
^ ^ lOOVK"^ ^ 4 much less than the system size N = 101. 
The displayed regions of high amplitudes show the positions 
where the particle localizes, (c) In the presence of dynamic 
disorder = 5 the quantum interference effects vanish. The 
main difference between dynamic disorder (c) and static dis- 
order (b) is that in (c) the regions with high values of P(x,t) 
keep changing leading to randomization so that the particle 
can still move but in a "classical" fashion. 



Anderson transition from extended to localized states is 
expected, via an intermediate chaotic regime which is 
rather better for quantum propagation. The probability 
density of Fig. 1(b) is shown to stay around the mid- 
dle site where the initial wave-packet has maximum am- 
plitude and it remains there for longer times. For the 
uniform initial state of Fig. 2(b) the larger amplitudes 
remain on many sites indefinitely. 



dynamic disorder 

A quantum walk can operate in the presence of dy- 
namic disorder but only for short times since for longer 
times its motion becomes entirely classical, indistinguish- 
able from an ordinary random walk. The effect of dy- 
namic disorder is equivalent to introducing coin chaos 
which makes the quantum coherence disappear [22]. In 
order to see this decoherence effect we have chosen a ran- 
dom eix^t) rapidly varying with t by: 

(i) Updating at random the diagonal site energies 
e(x,t) at a time length comparable and often much 
smaller to the time step of the numerical method. 

(ii) Varying the diagonal energies by 

e(x, t) = amp * cos{u;xt + (px) (3) 

where amp, u and (j) are the amplitude, frequency and 
phase for the motion of levels, respectively. We chose 
to vary the frequency oo at random uniformly within the 
interval [0, 27r] and fixed the phase (j)x to zero. 

From Figs. 1(c), (d) and 2(c), we can see the effect 
of dynamic disorder for the two initial wave-packets, S- 
function and broad, respectively. In Fig. 1(c) the quan- 
tum motion seen on the left hand side of the figure quickly 
disappears and this also happens in Fig. 1(d) where clas- 
sical diffusive motion is seen to arise. In Fig. 2(c) the 
randomization effects of dynamic disorder become more 
obvious and could be contrasted with quantum localiza- 
tion for static disorder (Fig. 2(b)). The quantum-to- 
classical crossover takes place after a characteristic time 
tc (X W-^. 



III. QUANTUM TO CLASSICAL CROSSOVER 
FOR DYNAMIC DISORDER 

decoherence in an A^-site chain 

The decoherence effect of dynamical disorder which 
turns the quantum wave propagation into classical dif- 
fusion is shown in Fig. 3 for a self-expanding chain. 
The probability density P(x, t) gradually changes from a 
shape displaying two ballistic peaks of the quantum wave 
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FIG. 3: The probability density is shown to display the grad- 
ual decoherence, from ballistic for small t to diffusive evolu- 
tion for long with strength of dynamic disorder W — 2. The 
two ballistic fronts in the quantum case gradually disappear 
and the shape approaches a Gaussian with classical diffusion. 



for small t towards a Gaussian for large t. In Fig. 4 the 
quantum to classical crossover is shown for the mean- 
square- variance cr^(t) and the autocorrelation function 
or return probability C{t) = ^ P{0,t')dt' . Eventually, 
the classical asymptotic laws cr^(t) cx t and C{t) oc 
set in after an initial period of quantum ballistic motion 
where cr^(t) cx and C{t) (xt~^. In Fig. 5 the effect of 
sinusoidal dynamic disorder is considered with constant 
amplitude and randomly varying the phase. The results 
are similar to Fig. 4 with the approach to the classical 
limit even faster in this case. From Figs. 4 and 5, ex- 
cept for no difference between the two types of dynamic 
disorder, we find that the crossover region between the 
ballistic law for small times and the diffusive law for long 
times is smooth having a mixed quantum and classical 
character. 




FIG. 4: (a) Log-log plot of the mean-square -displacement or 
variance cr^(t) vs. time t for various values for the dynamic 
disorder W. The crossover from ballistic quantum motion 
to diffusive classical motion occurs at earlier times as W in- 
creases, (b) The autocorrelation function or "return to the 
origin" probability C{t) vs. t for = 20 displays classical 
diffusive behavior. 



decoherence in a qubit 



The effect of dynamic disorder on the quantum evo- 
lution is displayed in the linear and log plots of Fig. 6 
which show the snapshots of the evolving spatial wave 
for an initial (5-function in a self-expanding chain. The 
decoherence effect of dynamical disorder is seen from the 
rapid approach to a Gaussian shape. The complete phase 
diagram is summarized in Fig. 7 with tc vs the strength 
of dynamic disorder W which displays a wide crossover 
grey color region where the law tc cx is approxi- 

mately obeyed. 



We have examined in detail the quantum walk in a 
random two- level system {N = 2) which has recently 
attracted attention in the context of quantum infor- 
mation processing. The operation of qubit and logical 
gates in the presence of a noisy environment is impor- 
tant for understanding quantum computers. The usual 
noise for such two-level system is usually due to vari- 
ous sources while non-Gaussian randomness arises from 
hopping background charges for different statistically in- 
dependent fluctuators. 
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FIG. 5: (a) The cr^(t) vs. time t same as in Fig. 4(a) but for 
the sinusoidal dynamic disorder of Eq. (4). (b) The C(t)) vs. 
time t same as in Fig. 4(b) but for the sinusoidal dynamic 
disorder of Eq. (4). 



time-dependent disorder W-2Q 




The time-dependent Hamiltonian is 
7 e2(t) 



(4) 



with random diagonal terms defined by 

{ait)) = 0, (e,(t)e,(i')) = SSijSit - t'),i,j = 1,2, (5) 

where S = W'^/12 measures the disorder chosen from a 
box distribution within [— W/2,+W/2]. The averaged 
matrix elements of the density matrix p can be obtained 
from a decoupling suggested in[23[ 

^p'li -ip22 7(^21 - PT2) (6) 

^Pi2 = -^idj^ + 7(P22 - pTT)) (7) 
and if we define 



FIG. 6: The snapshots for dynamic disorder of strength 
W=20 in a self-expanding chain where the approach to a 
Gaussian is seen, (a) The probability density P(x,t) = 
|^(x,t)|^ as a function of space x at fixed times t = 5000, 
t = 20000, t = 60000 and t = 80000. (b) The same as in (a) 
but for the log of the amplitude. 



the corresponding equations become 

p = -27 J, R = -2SR, J = -25 J + 27p. (9) 
Their general solutions (with appropriate constants) 



J 



C4 



A A 



-27 



A_ 



e^-^ (10) 
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FIG. 7: The quantum to classical crossover for dynamic dis- 
order of strength W occurs at tc oc W~^. This is shown by the 
grey area which was estimated from two sets of points con- 
nected via lines for the quadratic ballistic law to stop (blue 
line) and the linear diffusive law to begin (red line), respec- 
tively. 



R 



-S ± _ 4^2 



(11) 



by choosing as initial state one of the two levels with 
p(0) = l/2,i?(0) = J(0) =0. 

Finally, the averaged off-diagonal matrix element of 
the density matrix is easily shown to be 



Pl2{t)=lJ{t) 



■ sm 



( - 6H)e 



-5t 



(12) 

for 6'^ < 47^. If S'^ > 47^ in Eq. (12) the quantity un- 
der the square root changes sign and sin is replaced by 
sinh. Thus, for t ^ 00 the averaged density matrix ap- 
proaches half the unit matrix with only diagonal matrix 
elements and the quantum coherences described by pi2 
becoming zero, oscillating for the quantum case 6 < 2j 
and monotonically for the classical case (5 > 27. 

In order to examine the dephasing effect of dynamic 
disorder we have plotted in Fig. 8 the phase 6 of pi,2 vs. 
t obtained from numerical computations. Our results are 
presented for 7 = 1 and different values of disorder W 
which verify the critical value Wc = a/24 of the previ- 
ous analysis based on averages. The quantum coherence 
remains for weak dynamic disorder W < Wc while for 
higher dynamic disorder W > Wc the phase randomizes 
and the system becomes classical. 



IV. DISCUSSION 



Quantum walks are quantum analogues of classical 
random walks which have been proposed for quantum 
computation purposes to create quantum algorithms 
which run faster in quantum computers. They can also 
arise from mapping various physical problems (e.g. see 



FIG. 8: The phase of the off-diagonal matrix element of 
the density matrix pi, 2 vs. time for a two- level system with 
dynamic disorder = 0.1, 1 and 10. The decoherence due to 
dephasing is seen for the highest value W = 10 > Wc where 
completely randomizes. 



|24j). Some quantum algorithms which speed-up clas- 
sical methods have already been successfully employed, 
such as for search problems on graphs. These algorithms 
which show amplitude amplification during the evolu- 
tion could be efficiently implemented in a quantum com- 
puter. However, since they are often confronted with 
disorder we have examined how quantum wave-packets 
move in the presence of disorder, by computing the prob- 
ability density P{x,t) = |^(x,t)p from solving the time- 
dependent Schrodinger equation in the discrete space x 
of one-dimensional lattice. The static disorder due to im- 
perfections and the dynamic disorder due to the environ- 
ment could become obvious in scattering from nanostruc- 
tures or can appear from environmental noise, averaging 
over measurements, etc. 

Our main conclusions from the quantum evolution of S- 
function and very broad initial spatial wave-packets are: 
(1) In random media quantum walks can perform even 
worse than their classical counterparts since Anderson lo- 
calization from negative quantum interference completely 
stops the quantum walk although the corresponding or- 
dinary random walks in the presence of disorder can still 
move despite infinitely slowly. (2) For dynamic disorder 
we have no benefit from quantum walks either, since for 
longer times the ballistic evolution for small-t crosses over 
to classical diffusion for long-t and the quantum walks be- 
come classical via a quantum to classical crossover. (3) 
The answer to the question "what slows down the quan- 
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turn walk?" is, on one hand, "static disorder via nega- 
tive quantum interference" and, on the other hand, "dy- 
namic disorder at long enough times which slows down 
the quantum walk and makes it no different from or- 
dinary random walk". Therefore, quantum interference 
in random media can hold surprises for quantum walks 
and their advantages should appear for weak disorder or 
short times only. In higher dimensions quantum walks 
are also expected to operate for weak disorder to avoid 
quantum localization. In conclusion, our computations 
show Anderson localization or decoherence as the main 
enemies of quantum walks in the presence of static and 
dynamic disorder, respectively, which destroy their well- 
known quadratic or exponential speed-up. Our study 
could be useful towards creating better quantum search 
algorithms in the presence of disorder [25|. 
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